Entropy Production at RHIC 
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For central heavy ion collisions at the RHIC energy, the entropy per unit rapidity dS/dy at 
freeze-out is extracted with minimal model dependence from available experimental measurements 
of particle yields, spectra, and source sizes estimated from two-particle interferometry. The extracted 
entropy rapidity density is consistent with lattice gauge theory results for a thermalized quark-gluon 
plasma with an energy density estimated from transverse energy production at RHIC. 
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Lattice quantum chromodynamics (QCD) calculations 
indicate [l| that at zero net baryon density at a criti- 
cal temperature of T c ~ 170 MeV and energy density 
e c ~ 1 GeV/fm 3 a color deconfined and a chirally restored 
quark-gluon plasma (QGP) phase is formed. These en- 
ergy densities are attained in relativistic heavy ion colli- 
sions, where it is believed that a dense system of quarks 
and gluons is created which undergoes rapid collective 
expansion before the partons hadronize and eventually 
decouple. 

The ongoing experimental program at RHIC is de- 
voted to the search for the signals of QGP. Unfortunately, 
measurements are confined to final-state particles which 
are mostly hadrons. One interesting quantity that may 
provide valuable insight into the state of the matter in 
the early stages of the collision is the net entropy which 
is roughly conserved between initial thermalization and 
freeze-out 2, 3J. After freeze-out, when particles are free- 
streaming, the entropy is fixed. Entropy would also be 
conserved during the expansion stage if the mean free 
paths were very short. Viscous effects, shock waves, and 
the decoupling process can only increase the entropy. 
Thus, the entropy from the final state provides an up- 
per bound for the entropy of the initial state. 

Compared to a pion gas at the same energy density, 
the QGP is a high-entropy state as can be understood 
by considering the equation, 
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where AS is the change in entropy resulting from the 
addition of an amount of heat Q into a fixed volume 
at temperature T, Lattice calculations show that the 
transition is nearly a first order with a distinct peak in 
the specific heat near T c . In this region, the temperature 
stays relatively fixed while the energy density changes 
substantially. Thus, the existence of a latent heat keeps 
the temperature lower than a pion gas at the same energy 
density, which results in a higher entropy as can be seen 
from Eq. JTJ. At high energy densities, the temperature 
of a QGP remains lower than that of the pion gas due 
to the disparity in the number of degrees of freedom. 



Asymptotically, the energy density of a QGP approaches 
that of a massless gas of quarks and hadrons, 



e= N, 



(2) 



where Nb and Np arc the number of Bosonic and 
Fermionic degrees of freedom. Whereas a pion gas has 
3 bosonic degrees of freedom, a QGP has 47.5 effective 
degrees of freedom if one considers gluons and the three 
light quark flavors. Thus, relative to a pion gas, the 
QGP has a lower temperature and higher entropy. If the 
partonic degrees of freedom were not fully liberated at 
RHIC, one would thus expect a lower entropy for the 
same energy density. Furthermore, if one were to con- 
sider a full hadron gas, incorporating all the resonances 
from the particle data book at energy densities far above 
the point where the hadronic models are justified, one 
would expect to have a state of matter with even more 
entropy than a QGP for temperatures near or above 250 
MeV. 

The entropy density per unity rapidity dS/dy can be 
calculated from the final phase space density, /(p, r, tj), 
evaluated at the freeze-out time tf. Spectra provide in- 
formation only about the final momenta, p. The spa- 
tial information in / must be determined by correla- 
tions which can provide the characteristic dimensions of 
/(p,r,t) for given values of p 0,0- In fact, it has been 
suggested that the softer equation of state associated 
with a first order phase transition to QGP with a large 
latent heat results in long lived and/or large source sizes 
El2l0 0,03 The fact that the QGP models tend to ac- 
curately reproduce spectra while over-predicting source 
sizes H \m IH IH 03 suggests that the entropy at 
RHIC might fall below the expectations for a quark-gluon 
plasma. 

In this letter we extract dS/dy for hadrons at freeze-out 
from yields and spectra combined with the source dimen- 
sions in central Au+Au collisions at the RHIC energy 
of y/s = 130 A GcV. By comparing this estimated net 
hadron dS/dy with the entropy at the early stage in the 
parton medium, we address several key issues concerning 
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the initial stage of the collisions and the equilibration of 
partonic degrees of freedom. 

In the Landau quasi-particle approximation the en- 
tropy is given by 

S = (2J+1) |^^[-/ln/±(l±/)ln(l±/)] 

« (2 J + 1) / ^ [-/ In f + f± / 72] , (3) 

where the ± sign corresponds to bosons/fermions with 
total spin J. In the classical limit / « 1, the f 2 and 
higher orders of / can be neglected. The classical limit is 
sufficient for all particles except for pions. Even for pions, 
the f 2 term contributes only on the 3% level. Assuming a 
three-dimensional Gaussian form for the density profile in 
coordinate space, the single-particlephase space density 
at any time t at or after freeze-out ja, [w| is 



f{r,p,t) = / max (p)exp 
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where R x , R y and R z are implicit functions of p. Since 
spectra are obtained from integrating the phase space 
density, /max, and therefore /(p, r, t), can be determined 
from the radii and spectra. 



dN 
d 3 p 

/max(p) 



(2J+1) /(^ /(r ' P) 

(2tt) 3 / 2 dN 1 
(2J+l)d 3 p'R 3 ' 



(5) 



where R 3 is the product of the three radii. Chemical 
equilibration has not been invoked as the expressions do 
not involve temperature or chemical potentials. Instead 
the phase space densities and therefore the entropies will 
be functions of the spectra and the source sizes which 
can be inferred from experiment. 

In Hanbury-Brown Twiss (HBT) correlation measure- 
ments for identical particles, if the single-particle emis- 
sion source is characterized by a Gaussian as in Eq. 
the Gaussian widths are each inversely proportional to 
the source dimensions in the canonically conjugate spa- 
tial variables 0] . One then obtains the two-particle cor- 
relation function in terms of the radii, 



C(K,q) = 



P( P l,p 2 ) 
P(P1)P(P2) 

1 + exp (-q 2 Rl - q 2 R 2 - q 2 R 2 ) , (6) 



where K = (pi+p2)/2 and q = pi — P2- The momentum 
difference q has been decomposed into "out-side-long" 
components defined by the beam direction qi, parallel to 
the total momentum of the pair in the transverse plane 
q 0l and the direction orthogonal to the above two direc- 
tion q s 0, 0- Data are generally analyzed in the lon- 
gitudinal co-moving system (LCMS) frame in which the 
longitudinal component of the pair momentum vanishes. 



The Gaussian widths of Eq. © in the LCMS frame 
are related to the radii in the two-particle rest frame via 
a contraction factor 7 = Ex/m. Phase space densities 
depend on the product of the three radii. In the rest 
frame of the pair, this product is referred to as Rf nv = 
jR R s Rl. Substituting the HBT radii in Eq. ©, the 
entropy rapidity density from Eq. @ in the classical 
limit is then 
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Thus, dS/dy of the final state can be estimated from the 
Gaussian source-size parameters and spectra. 

The two-pion correlation function has been measured 
by STAR 16] and PHENIX [HI for ~ 11% central 
Au+Au collisions at y/s = 130A GeV at midrapid- 
ity In spite of about 70% increase in pion multiplic- 
ity from SPS to the RHIC energy the extracted radius, 
i?i nv ~ 6.3 fm at RHIC, was very similar to what was 
measured at the SPS 18]. Moreover, the prediction 
of R larger than R s as a QGP signature was not ob- 
served. Instead R /R s w 1 was measured over a broad 
range of transverse momenta Kt = (pTi + Pt 2 )/2- The 
smaller values of Ri and R , as compared to hydrody- 
namic models, may be suggestive of a rather high final 
phase space density I^Lllfll]. and therefore, a low entropy. 
In our estimate of 7r~ entropy we use the scaling rela- 
tion Ri = Aij \/rriT where mj- = yj K T + m 2 with the 
fitted value A t = 3.32 ± 0.03 fm GeV 1 / 2 17 : ). While for 
the s-o direction we employ the blast-wave inspired pa- 
rameterization, R 2 = R 2 = i?g 0om /(l + /3 2 tot/T), where 
■Rgeom = 6.7±0.2 fm, j3 = 0.69 and the freeze-out temper- 
ature is T = 125 MeV ^3]- The measured corresponding 
single particle ir~ spectra at RHIC (after correcting 
for electromagnetic and weak decays) is fitted to a Bose- 
Einstein distribution EdN/d 3 p = A' /[exp(m T /T oS ) - 1] 
where A' = 739 GeV -2 and effective (slope) temperature 
T Q g = 197 MeV. The entropy per unit rapidity estimated 
from Eq. Q is dS/dy\^- = 680, and is shown in Table 
I. 

Pions resulting from electromagnetic decays, i.e., those 
from 77 and rj' decays, should be subtracted from the spec- 
tra. From a thermal model with a chemical equilibrium 
temperature of 170 MeV, we estimate that 12% of the 
pions come from such decays and correct the spectra ac- 
cordingly. We add in the entropy of the rj and rj' by 
assuming that the 77 contributes the same amount of en- 
tropy as a single species of kaon which has a similar mass, 
and that the 77' contributes half the amount of entropy 
as the average of protons and anti- protons. The factor of 
1/2 is applied to account for the lower spin degeneracy 
of the vl . The entropies of the r\ and r( are displayed in 
Table I. 
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TABLE I: The entropy content per unit rapidity, dS/dy, for hadrons produced at midrapidity in 
collisions at the RHIC energy of yfs — 130A GeV. 



11% most central Au+Au 



7T 


v(v') 


K~ (K+) 


P (p) 


A+£° (A+E°) 


E* (E±) 


& (s ± ) 


Q. + (n-) 


680 


256 (69) 


256 (286) 


118 (159) 


81 (111) 


31 (42) 


23 (27) 


5 (5) 



For a dilute gas in chemical equilibrium, it was demon- 
strated 0] that the entropy per nucleon carried by the 
nucleons is related to the ratio of the deuteron to the 
proton yield, 



S = 5/2-lnR dp /3V2 



(8) 



If, within a given phase space cell, the ratio of a species 
a to the product of two species b and c, which can un- 
dergo a reaction a + b <-> c, is determined by the ratio 
of available phase space for bound and unbound states, 
one can determine the Gaussian source parameters from 
coalescence arguments |2l|. 



R\, 



r 3/2_ 



(2 J c + 1) 



(2J a + l)(2J b + l) 
E a dN a \ ( E b dN b 



E B dN, 



TO d 3 p a J \rrib d 3 pb J I \ m c d 3 p 



(9) 



This expression is based on the assumption that the rel- 
ative populations within a given phase space cell are de- 
termined by the ratios of the available phase space for 
bound and unbound states; whereas Eq. (JSJ requires the 
assumption that all phase space cells are equilibrated ac- 
cording to a single chemical potential. Once -R; lw is ob- 
tained as a function of the momentum, one can extract 
the entropy from Eq. (JJJ. 

Assuming that p and n possess identical phase space 
densities, which in turn is identical to that for d at the 
same relative velocity, one can determine a source size 
for p using measured anti-proton and anti-deuteron spec- 
tra. Figure n shows the R mv for antiprotons (open trian- 
gles) estimated (using Eq. (JSJ) from the measured (2^ 
spectra of d and weak-decay-corrected ps at three Kt 
bins for the ~ 18% most central Au+Au collisions at 
■s/s = 130A GeV. As expected, the radii for the massive 
p exhibits a smaller decrease with Kt compared to that of 
pions. The p transverse momentum spectra for this cen- 
trality from the STAR ^ and the PHENIX [U data are 
corrected for the weak hyperon decay [24j and matched 
in the overlap region. We fit this spectra in the range 
0.25 < px < 3.2 GeV/c with a Gaussian distribution 
EdN/d 3 p = dN/dyexp(-p T /2a 2 )/(2TTa 2 ) resulting in a 
yield of dN/dy = 11.4 and width a = 872 MeV. Given 
the limited range and poor statistics for anti-deuterons, 
we choose a constant average value, Ri nv = 5.34 fm, in- 
dependent of Kt- Applying Eq. (JJJ), the entropy carried 
by the anti-protons is then dS/dy\p = 118. 



The procedure to estimate the p source radii can be 
identically applied to extract the kaon radii using the <j> 
meson, «-> K + K~, in place of the deuteron. The co- 
alescence arguments used in Eq. @ ignore the binding 
energy of the composite particle, which is justified if the 
binding energies are much less than the breakup temper- 
atures. This is an excellent assumption for the deuteron 
and a reasonable assumption for the <f> which is a reso- 
nance with a binding energy of B = — 30 MeV. Although 
we neglect the binding energy in this analysis, one could 
adjust the expression for Rf„„ by a factor of e B / T to ac- 
count for the binding energy |21| . 
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FIG. 1: Top panel: The K T dependence of HBT radii for 
K~ (open circles) and p (open triangles) estimated from in- 
variant yields at midrapidity for central Au+Au collisions at 
y/s — 130A GeV. The measured ir~ radii (solid circles) from 
PHENIX [Tt| are also shown. Bottom panel: The average 
freeze-out phase space density as a function of Kt for the 
hadrons. 

Figure^displays the K~ invariant radii extracted from 
measured midrapidity cj> [25| and K~ j2(| yields for the 
11% most central Au+Au collisions at y/s = 130 A GeV. 
As expected the kaon radii fall between those for the pro- 
ton and pion. The Kt dependence for K~(K + ) radii is 
parameterized as i?; nv = A/^/my, with A = 4.45(4.74) 
fm GeV 1 / 2 . The estimated kaon radii are consistent 
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with that of the preliminary three-dimensional measure- 
ments 27] at the highest RHIC energy of ^/s = 200A 
GeV. We note that the ttjt dependence of the HBT radii 
i?i„ v = A I ^JrriT for these hadrons at RHIC are not con- 
sistent with one scaling coefficient as observed at the SPS 
energy [Hj]; the fitted values here are A = 3.81, 4.45, 5.32 
fm GeV 1 / 2 for ir~,K~,p. The measured single par- 
ticle nvr spectra for K~ and K + [2(| are fitted with 
an exponential function EdN/d 3 p = dN / dy exp[— (hit — 
TrbK)/T]/(2-KT{rriK + T)) that results in a rapidity den- 
sity of dN/dy = 38.2(42.1) for R-{K+) using a slope 
parameter of T = 280 MeV. Inserting the spectra along 
with the the parameterization of i?i nv into Eq. J7J) yields 
dS/dy\ K - — 256 and dS/dy\ K + = 286. We assume that 
the neutral (anti)kaons carry the same entropy as the 
corresponding charged (anti)kaons. 

The measured [24| spectra for A for the 11% most 
central collisions at midrapidity has feed-down from the 
electromagnetic decay of the £ and the weak decays of 
the S + , S , and £1 + . The feed-down from the charged 
S and Q states is corrected by using the exponential 
parametrization of the the measured spectra for S and 
il + of Ref. . The corrected A spectra is then fitted 
with an exponential function resulting in dN/dy = 7.4 
and slope parameter T — 390 MeV. The entropy rapid- 
ity density for the antihyperons estimated by assuming 
that they have the same source radii as the antipro- 
tons are presented in Table I. To estimate dS/dy for 

£ , we note that the particle density in the statisti- 
cal model in the Boltzmann approximation is N/V — 
gm 2 TK2(m/T) exp(/i/T) / (2tt 2 ), where the notations are 
as usual. Assuming isospin symmetry E + « E° w S , 
and same temperature « Tg- = 173 MeV, the ther- 

modynamical model gives the ratio A/S = 1.47 and 
thus = (A + E°)/2.94. The dS/dy extracted under 
this assumption for £ is listed in Table I. 

To extract the entropy density for baryons, we employ 
the measured ratio for the yields p/p = 0.72, A/A = 0.71, 
= 0.83, H + /rr = 0.95 [13, and assume that all 
(anti)baryons have the same source radii as antiprotons. 
Furthermore, the ratios are found to be independent over 
a large range a pr suggesting that the transverse momen- 
tum distribution and thus the slope temperatures of the 
baryons and their antibaryons can be taken to be identi- 
cal. The dS/dy for the baryons thus estimated are shown 
in Table I. 

The values of i?i nv shown in Fig. ^ (top panel) may 
be used in conjunction with the corresponding spectra of 
the hadrons to calculate at a given px the average phase 
space density at freeze-out (J(pt)) — ^ r /2v / 2, where F 
is given by Eq. (JJJ) . This is displayed in Fig. ^ (bottom 
panel) where the pion (/) at RHIC is remarkably higher 
than the values extracted at SPS [13, |3jJ and AGS |32| 
energies. 



Under the assumption of isospin symmetry the total 
entropy rapidity density for the hadron-gas at freeze-out 
for central collisions at RHIC is 4451. Despite the large 
uncertainty in source sizes, this number is surprisingly 
robust since the radii only enter logarithmically. A 30% 
change of Rf nv changes the entropy by 0.26 units per par- 
ticle. For baryons, which have approximately 10 units of 
entropy per particle, the overall entropy would change 
by only a few percent. Another source of error comes 
from estimating the contribution of the rj and rf. How- 
ever, since an r\ carries two to three times more entropy 
than a pion, and since the rj decays into two or three 
pions, the net entropy is rather insensitive to the frac- 
tion of pions assigned to 77s. The same would be true if 
one were to explicitly include the entropy of the <j> meson 
while correspondingly reducing the kaon yields. Since the 
yields are not measured to better than 10% accuracy, the 
systematic error in estimating the final-state entropy is 
estimated ~ 10%. 

We now estimate the entropy per rapidity density gen- 
erated in the QCD partonic matter in the early stage at 
RHIC. Assuming boost-invariant longitudinal expansion, 
the energy density averaged over transverse coordinate is 

M 



eBj(r) = 



dEi 



1 



dy TirR 2 



(10) 



Using dE T /dy » 542 GeV |3J and wR 2 = 120 fm 2 for 
the ~ 11% most central Au+Au collisions, one obtains 
e(r) = 4.5 GcV/fm 3 at a proper time r = 1 fm/c. Present 
data from lattice QCD simulation at /is = shows a 
rapid increase of e/T by about an order of magnitude 
within a narrow temperature range AT < 0.1T C ~ 16 
MeV near the critical temperature T c Q . This may indi- 
cate a first order transition, or nearly first-order transi- 
tion. For three degenerate quark flavors, rif = 3, with a 
mass m/T — 0.4, transition to a plasma phase occurs at 
T c = 154 MeV at a critical energy density e/T c 4 ~ 6. As- 
suming the initial matter at RHIC is already thcrmalizcd 
at 1 fm/c, lattice data indicates that e(r) = 4.5 GcV/fm 3 
would correspond to an (average) initial temperature of 
T « 230 MeV and pressure density p w 3T 4 = 1.1 
GeV/fm 3 . The resulting entropy density for the 3 fla- 
vor system is s = (e + p)/T = 24.3 fm~ 3 , and the average 
entropy rapidity density at r = 1 fm/c is 



<iS_ 

dy 



QGP;n /= 3 



d 2 X T S ~ 2922 



(11) 



The Bjorken estimate of the energy density only pro- 
vides a lower bound since it ignores the energy of the 
local longitudinal (along the beam axis) motion and ne- 
glects the longitudinal work caused by the expansion of 
the plasma. The inclusion of local longitudinal motion 
would increase the estimate by a factor of 4/ir if the 
emitted particles were massless and if there was no trans- 
verse collective flow. Accounting for transverse flow and 
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FIG. 2: The entropy per unit rapidity from lattice calculations 
(solid line) is displayed as a function of the eo, the energy den- 
sity at r = 1 fm/c. The horizontal band shows the final-state 
entropy extracted from experiment. Hydrodynamic models 
typically have energy densities nearing 10 GeV/fm 3 , which 
would push the limit for entropy from the final state. 



masses should result into the increase being in the 10- 
15% range. For massless particles undergoing a purely 
one-dimensional Bjorken expansion, the energy density 
would fall as r -4 / 3 , and dEr/dy would change by a fac- 
tor of r -1 / 3 . This would be a factor of 2 for times be- 
tween 1 and 8 fm/c. The loss of Et from longitudinal 
work is proportional to the irR 2 J Pdr in hydrodynamic 
models, or irR 2 T zz dr if the stress tensor is non-isotropic, 
i.e., viscous effects are included. The effects of viscosity 
have been studied within the context of cascade models. 
Whereas, onc-dimcnsional hydrodynamic models might 
show that dEx / 'dy is higher than that at breakup, the 
inclusion of viscosity reduces that increase to the 50% 
range |35j, |36| . After including the effects of longitudinal 
work and longitudinal motion, estimates of the initial 
average energy density might be in the range of 7 ± 2 
GeV/fm 3 , depending on the equation of state and the 
viscosity. 

Of course, the initial entropy is less than the final 
state entropy. One expects the entropy to increase dur- 
ing the reaction by the order of 10% from viscous effects 
or shocks 37] . If one were to lower the final-state esti- 
mate of dS/dy by 10% it would come close to the value of 
dS/dy predicted by the lattice gauge calculation. A lack 
of complete thcrmalization of the initial state would also 
point to a lack of entropy at early times. If the initial 
energy was stored in classical gluon fields j^, the en- 
tropy would be negligible. An undcrpopulation of quark 
degrees of freedom would also lower the entropy. Since 
flavor and color charges cannot sample the entire volume, 
finite-size effects can also reduce the initial entropy, but 
only by a few percent |3flt |40j . Since the amount of en- 
tropy produced during the expansion is unknown, it is 



impossible to rule out that the initial state might have 
been of lower entropy. 

Thus, our estimate of the final-state entropy, dS/dy = 
4451, does not seem inconsistent with expectations from 
lattice calculations from analyzing Fig. It therefore 
seems puzzling that hydrodynamic calculations based on 
lattice-inspired equations of state matched spectra and 
failed to fit HBT sizes H E3, GH E3 An overestimate 
of a factor of two in Rf nv would suggest an overestimate 
of the entropy of the order of ln(2) per pion, or ~ 400 
units per rapidity. However, it appears from the data 
that other degrees of freedom, especially the baryonic 
degrees of freedom, compensate for the lack of entropy 
in the pions. Although the pions comprise two thirds 
of the final-state particles, they provide less than half 
the final-state entropy. The shift of entropy from pions 
to the (anti)baryon sector was discussed by Rapp and 
Shuryak [4l| as a natural consequence of an increased 
phase space density of pions. This follows from the reac- 
tion pp <-* Ntt. Since pion number population is difficult 
to maintain, the overpopulation of pions was expected 
in an iscntropic expansion [42j. An additional factor in 
increasing the overpopulation of pions would be a low- 
ering of hadronic masses associated with chiral symme- 
try restoration 43] . A lowering of hadron masses would 
lead to an increase in the population of heavier hadrons, 
many of which would decay back to multiple pions after 
the normal QCD vacuum is restored |44| . 

It should be stressed that the small pionic sources ex- 
tracted from HBT remain a puzzle. Since the overall en- 
tropy of the final state is consistent with the lattice gauge 
theory, the possibility remains open that a quark-gluon 
plasma is indeed created at RHIC but that the dynam- 
ics of the expansion and decoupling somehow differ from 
current hydrodynamic descriptions [iH |46| . 

To illustrate the sensitivity of the entropy to the num- 
ber of degrees of freedom, Fig. [21 also displays dS/dy for 
a pion gas. Since the pion gas has only three effective 
degrees of freedom, it has a much higher temperature 
for a given energy density than does the lattice equa- 
tion of state. From Eq. it therefore has a lower en- 
tropy when compared at the same energy density. At 
high temperatures, a hadron gas would effectively incor- 
porate more degrees of freedom, and for energy densities 
near or above 10 GeV/fm 3 would have more effective 
degrees of freedom than a plasma due to the extremely 
large number of baryonic resonances with masses between 
1.0 and 2.0 GeV. However, since the volume of a baryon 
is in the neighborhood of 2 fm 3 , energy densities above 
0.5 GeV/fm 3 preclude a serious description in terms of a 
non-interacting hadron gas. 

In conclusion, we have calculated the final-state en- 
tropy from 130^4 GeV Au+Au collisions at RHIC to be 
dS/dy — 4451. We estimate that this number, which is 
based solely on measured spectra and extracted source- 
size radii of the principal final-state constituents, to have 
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a systematic uncertainty of 10%. This number is not out 
of line with expectations from lattice gauge theory. Im- 
proved experimental measurements of spectra and source 
size might reduce the systematic uncertainties to the 5% 
level, but it is the uncertainty in the reaction dynam- 
ics between the initial and final stages of the collision 
that precludes our making a more definitive conclusion 
about the initial state entropy. Although our number 
represents an important upper bound on the initial state 
entropy, the amount of entropy generated throughout the 
reaction and the initial energy density remain objects of 
conjecture. These issues will be settled largely by other 
observables. For instance, elliptic flow measurements are 
sensitive to the viscosity of the initial stages 0] and 
HBT measurements provide important three-dimensional 
constraints about the dynamics. 

This work was supported by the National Science 
Foundation under Grant No. PHY-0070818, and by the 
U.S. Dept. of Energy, Grant DE-FG02-03ER41259. 
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